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ABSTRACT 


Tim general characteristics of pulsar Y-ray spectra are presented for 
a model where the Y-rays are produced by curvature radiation from energetic 
particles above the polar cap and attenuated by pair production. The shape 
of the spectrum is found to depend on pulsar period, magnetic field strength, 
and primary particle energy. By a comparison of numerically calculated 
spectra with the observed spectra of the Crab and Vela pulsars, it is 
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determined that primary particles must be accelerated to energies of ~ 3x10 
o 

me . A general formula for pulsar Y-ray luminosity is determined and is 

found to depend on period and field strength. The Y-ray efficiency scales 

18 7 

with characteristic ago as T up to r * 3-4x10 yr», whore there may be a 
turnoff in the Y-ray emission. 
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I, INTRODUCTION 

High energy pulsed Y-ray emission from the two shortest period radio 
pulsars is now well established observatlonally. The Crab and Vela pulsars 
have been studied extensively in the last few years by SAS2 (Knlffen et al. 

1974, Thompson et al, 1975) and by COSB (Bennett et al, 1977) at energies 
above 30 MeV. The Y-ray emission from these pulsars io a much higher fraction 
of the total energy loss rate than the radio emission and, therefore, must 
play a major role in the radiation mechanism. The energy spectra of the 
Crab and Vela have recently been reanalyzed in the range 50 MeV-2 GeV from 
COSB data (Llchtl et al. 1980), and they fit steep power law*i with slopes of 

i * 

-2,17 and -1.89, 

Although only a small fraction of the 328 known radio pulsars havs 
been detected in the Y-ray range, it is possible that many of these pulsars 
are emitting pulsed Y-rays below the present detectability limits. If Y-ray 
emission is a general characteristic, an important question is whether the 
production mechanism is the same for all pulsars. The similar pulse shapes 
and spectra of the Crab and Vela argue strongly in favor of a single mechanism 
for these two pulsars, though their properties are quite different at other 
energies. From a theoretical point of view, it is essential to try to establish 

the characteristics of pulsar Y-ray emission in specific models for the whole 

* 

range of pulsar parameters such as period and magnetic field strength. When 
Y-ray data on other pulsars becomes available from more sensitive detectors, 
it may then be possible to use these characteristics as a test of certain 
models for Y-ray pulsars as a class. 
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Detection of pulsed emission from PSR1747— 46 has also been reported 
(Ogelman et al. 1976, Thompson et al. 1976), but this result is now in question 
because of a recent radio timing observation (Newton et al„ 1980) which 
measured a much lower period derivative than was used to determine the Y-ray 
light curve . 


It is the purpose of this paper to outline some observable charac- 
teristics of pulsar Y-ray emission In a particular model and to determine 
a general luminosity formula which is a function of pulsar parameters. 

This modal assumes that the Y-ray emission above 100 MeV is produced by 
primary particles undergoing curvature radiation in the magnetic field 
above the pulsar polar cap. These particles have been accelerated along 
field linos by the intense electric fields induced by rotation. Because 
of the extremely high energies to which the particles are accelerated 
(Sturrock 1071; Ruderman and Sutherland 1075; Scharleraann, Arons and Fawloy 
1978), their curvature radiation spectrum extends well into the Y-ray range. 
The most important attenuation mechanism affecting these Y-rays is magnetic 
field induced pair production, which becomes important around 500 MeV. 

The details of this model were presented in an earlier paper (Harding, 
Tademaru and Esposito 1978, hereafter referred to as Paper X), in which 
the Y-rays were traced from their emission points through the rotating 
magnetosphere and attenuated by pair production. Optical depths, pulse 
shapes, and spectra were calculated in the case where the emitting 
primary particles are monoenergetic. It was found that pulsar rotation 
increases the amount of pair production, primarily through the induced 
electric field perpendicular to the magnetic field. The electric field 
increases the pair production rate above the value it would have in a 
pure magnetic field (Daugherty and Lerche 1975). As a result, high 
energy Y-rays are more severely attenuated in short period pulsars. Salvati 
and Massaro (1978) obtained similar results using the same type of model, but 
including some synchrotron radiation from secondary pairs. 
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This paper assumes a monoenergetic Injection spectrum of primaries and 
takes into account the energy loss of the particles as they radiate (S II), 
This produces steeper Y-ray spectra than, those calculated in Paper I for 
constant primary particle energy. Numerically calculated spectra for various 
pulsar parameters and primary particle energies are presented in § III and 
are compared to the COSH spectra of the Crab and Vela. In S IV, the 
dependence of Y-ray luminosity on period and magnetic field strength is 
derived from the integrated spectra. A prediction of the dependence of 
efficiency on characteristic age is also discussed. 

II. CURVATURE RADIATION SPECTRUM 

paper I describes the basic assumptions and approximations of the 
model and also the numerical technique used to calculate the pair production 
rate as the photons move out from their emission points. The emission 
spectrum here differs from that described in Paper I in that it allows for 
energy loss by the emitting particles as they move out from the accelerating 
region. The calculation does not deal with the specifics of the acceleration 
mechanism but rather assumes that the particles somehow reach highly 
relativistic energies near the star. 

Several methods of particle acceleration above the polar cap have 
been discussed (Ruderman and Sutherland 1975; Scharlemann Arons and Fawley 
1978; Arons and Scharlemann 1979) as well as acceleration in the outer 
magnetosphere (Cheng, Ruderman, and Sutherland 1976). These studies determine 
the energy, Y Q » to which the particles are accelerated as a function of 
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period P, magnetic field strength li Q » and radius of curvatuce, R c , of the 

field lines i It is Pound that, while very large electric fields can be induced 

along the B field lines, acceleration is basically limited by pair production 

discharges which short out the E field. As a result of these competing 

processes each of which depend on P and B Q , the net potential drop has only a 

weak dependence on these parameters, In this calculation, we will simply 

assume that Y„ does not depend either on P or B . 

o o 


We assume that a monoenergetic beam of primary particles of energy Y 0 
is injected uniformly over the polar cap and moves out along dipole field 
lines. As the particles radiate, they will be losing energy at a rate. 
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where R c j | is the magnetic field line radius of curvature for a dipole 
field. The equilibrium distribution of particles at (r,0) will be, 


r N o ( c /*c ) ,*o* 

N(Y ’ r,0) W/dt) curv ( r > ’ 
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where Y max (r,0) is the maximum particle energy at (r,0), obtained by 
integrating the expression for (dY/dt) curv along a field line; 


V max (r ’ 0) 


Y q [3.16x 10" 13 | ln(~ ) Y* + U~ X/3 


(3) 
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N la the density of primary particles at the stellar surface, r , emerging 
o u 

from the acceleration zone. Since we are ignoring the details of the 
acceleration, we simply assume that N q is proportional to the corotation 
charge density: 
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where is the surface magnetic field in units of 10 gauss. The emission 
coefficient at point (r,0) will therefore be 


e(u), r, 0) - J* • N (Y, r, 9) V <u>) dY 
'max (r,8) 


(5) 
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(Jackson 1962) is the instantaneous power emitted per electron and U) c is 
the critical frequency, above which the spectrum cuts off. The upper 
integration limit in equation (5), Y f > is the lowest particle energy which 
can produce photons of frequency tu, or 


r 2 r 0 v /a 

<“> - it—" 1 / *■ 




Performing the Integration in equation (3), we find, 
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This is the energy spectrum radiated by primary particles in a volume 
element of the magnetosphere. The radiation will be highly beamed in a 
direction tangent to the field lines, and emission points are chosen so 
that photons in a given part of the pulse arrive at a stationary observer 
simultaneously. The determination of emission points and the subsequent 
attenuation of the radiation by pair production has been detailed in Paper 
X. Since the secondary pairs are produced at an angle to the field, they 
will lose most of their energy by synchrotron radiation. While this 
radiation may contribute to the Y-ray spectrum, the contribution above 
100 MeV will come from the high energy tail of the steep primary spectrum. 

We assume that this is a small effect, though it has not been fully determined, 
and neglect the radiation from secondary particles in the present calculation. 


III. CALCULATED SPECTRA 


Spectral points are calculated in this model by integrating the intensity 
at a given photon energy across the pulse profile. The free parameters 
determining the shape and magnitude of the spectrum are the period P, the 
surface magnetic field strength B q , the initial primary particle energy y , 
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and number density n, The angle between the magnetic field and spin axes is 
taken to be 90° in all calculations, with the observer's line of sight passing 
through the magnetic pole. 

For given values of P and B Q , Y 0 uniquely determines the shape of the 
spectrum and n determines the normalization. Figures 1, 2, and 3 show 
calculated spectra ter a - 1 and various values o £ B 0 , y Q and ?. The value 
of Y q clearly has a large effect on the curvature radiation cutoff which 
is in part responsible for the turnover in the spectrum at high energies. 

Pair production also steepens the high energy part of the spectrum and 
depends on both P and B Q , The effect of P can be seen by comparing the 
shapes of the spectra shown in each Figure, which are calculated for the 
same values of B Q and Y Q . The spectra steepen at lower energies for shorter 
period pulsars as a result of the higher pair production rates caused by 
rotation. Increasing B Q also increases pair production and, therefore, 
also has a steepening effect on the spectrum. The spectra for P ■ .089 s. 
in Figures 1 and 2 are nearly identical in shape to the spectra obtained by 
Salvati and Massaro (1978) in their Figures 4a and 4b for B q - 10 ,ff gauss. 

Since they have included a contribution from secondary pairs and our calcula- 
tions have not, it appears that the secondary pairs do not contribute 
significantly to the spectrum above 100 MeV. 

12 

If a value B Q * 10 gauss is assumed, then the best fit spectrum to 
the COSB data points for the Crab with P *• .033 s. is the one having 
Y q - 3xl0 7 and requires that n * 1.2 x 10 3 for normalization at 100 MeV (assuming 
a distance of 2 kpc.) A lower value of Y Q will fall below the high energy 
points. These same values for Y Q , n and B q , but with P - .089 s, also 
fit the spectral points for Vela quite well (assuming a distance of 0,5 kpc.) 
These spectra and their fits to the data are shown in Figure 4. The fact 
that both the Crab and Vela spectra can be fit with the same value of Y 

o 
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strongly supports the prediction of acceleration theories that the paifticle 
energy depends only weakly, if at all, on the pulsar period. It alsM indicates 
that, in order to account for the highest energy y-radiacion from these pulsars, 

the acceleration mechanism must he capable of producing particles with Y > 

7 3 ° 

3x10 , and particle densities n ~ 10 times the corotation charge density. Energies 

this high are not acheived in the present models for acceleration above the polar 

12 

cap where potentials only reach 10 volts, though particle densities of this 

* 1 

magnitude are. Possibly, the highest energy particles ere accelerated, or reac- 
celerated, in outer gaps where potentials can reach 10^ volts (Cheng, Ruderman 
and Sutherland 1976). As we have found in Paper X, high energy photons which are 
observed must actually be pr©# Jt'fid further out In the magnetosphere where the 

ootical deaths are lower, and the field line curvature is greater, 

* 

IV. Y-RAY LUMINOSITY 

The calculated spectra have been integrated for a range of pulsar periods 
and magnetic field strengths to obtain the dependence of Y-ray luminosity on 
these parameters. A fit to the numerically integrated spectra gives an 
expression for the luminosity above 100 MeV, normalized to the Crab with 
Ly (> 100 MeV) - 1.4 x 10 3 ® photon s“^ and ■ 3.8, 

Ly (> 100 MeV) - 1.2 x 10 35 B^ 5 P“ 1,7 photons s“l (8) 

IQ 

where B 12 is in units of 10 gauss. Tills formula gives Ly (> 100 MeV) - 

2.3 x 10 37 s-- 1 * for the parameters of Vela (B^ 2 " 3.4, P « .089 s), which 

37 

is in good agreement with the observed value 2.4 x 10 s"" 1 (Lichti et al. 1980). 
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where r Is the polar cap radius. If t Q is Independent of P and B Q , as we 
havo assumed, N 0 « B^' 1 , and r p - t 0 (O r o /c>* « ?"*, tbeal^ « B 0 P' 2 . 

The departures from this dependence In equation (8) therefore result from 
the inclusion of pair production attenuation. Since pair production Is 
greater In shorter period pulsars! their luminosities are decreased relative 


> 

to longer period pulsars. 




Luminosities and fluxes for other pulsars can be predicted from equation 
(8). The magnetic field strength Is usually estimated from the amount of 
radiation reaction torque implied by the value of the period derivative, P 
(Ostriker and Gunn 1969) 


B “ 3.2 * 10 19 &-■ - . ■/ < P 9)4 «*““• 
° 10“ g cm ^ 


(10) 


While It is reasonable to assume that all neutron stars have comparable 

masses, and therefore similar moments of inertia, I, and radii, r , It is 
r « 

also possible that the range of observed PP Indicates a range of neutron 

star masses. In other words, it is possible that the value of 3 q is the 

same for all pulsars. However, observational determinations of the masses 

j of neutron stars in orbit around binary companions (Rappaport and Joss 1979, 

i Taylor et. al. 1979) have all been consistent with 1.4 ± 0.2 M 0 , which seems 


j. 
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.§ 

to suggest a vary small maaa range, The observed values of PP also give 

10 13 

a reasonable range for B Q <2 x 10 - 2 x 10 gauaa) If a constant value 

g 

of I/r ” Is assumed, 
o 

a 

For pulsars with measured values of P, the estimates of B q , assuming 

45 2 6 

canonical values of I - 10 g cm and r Q « 10 cm, are used in deriving 
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the luminosity. For all other pulsars, a value of B q ■ 10 gauss is 
assumed. The pulsars with the highest predicted Y-ray fluxes are listed 
in Table 1 with tneir luminosities and other measured parameters from 
Manchester and Taylor (1980). Most of these have fluxes well below those 
of the Crab (PSR 0531 + 21) and Vela (PSR 0833-45), but a few have comparable 
predicted fluxes. Eight are above the point source detectability limit of 
COSB, but all are either outside the latitude range of the search or near 
the galactic center where the diffusa background is high. All of the pulsars 
in Table 1 with previously reported period derivatives have been searched for 
pulsed Y-ray emission by SAS2 and COSB. The predicted fluxes for PSR 0950 + 08, 
PSR 1642-03, PSR 1706-16 and PSR 1929 + 10 exceed the flux upper limits reported 
by SAS2 (Ogelman et al. 1976), but by no more than a factor of 3. This is not 
a serious discrepancy, considering that the errors in the SAS2 limits are of 
the same magnitude. 

Pulsar Y-ray efficiency is defined as the fraction of the total energy loss 
rate which is radiated in Y-rays. From equation (8) and the expression for the 
rotational energy loss rate, 

■* 

"if " * n2 1 ^3 erg 8-1 » (U) 

the efficiency above 100 MeV will be 

Ly (> 100 MeV) 2.0 x 10 31 ,95 p 1 * 3 

" ^“dE7dt 4^— Bl2 T~ 


( 12 ) 
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For constant moment of Inertia, the efficiency can be expressed In 

• i. 

terms of the pulsar's characteristic age, T - ^ (P/P)®, using equation (10): 

T) y - 4x10"^ V- 8 (13) 

* 15 1 

where P is In units of 10 s s~ A and t Is In years. For the two pulsars 

with measured Y-ray luminosities, there seems to be an increase in efficiency 

*— 1 3 

with age. Figure 5 shows that the quantity 1] P * plotted against t for 

these pulsars is consistent with a r^» 8 dependence. Equation (13) predicts 

that 7|y will exceed 1 above some characteristic age which is equal to ~ 3x10 7 
*' 15 1 

years for P ~ 10” A s s" . This suggests that the Y-ray emission mechanism 
may turn off before this point since the efficiency cannot exceed 1. One 
explanation for a Y-ray turnoff is that the primary particle energy, which 
was assumed to be constant in deriving equation (13), begins decreasing 
with increasing P at the point where the induced electric fields are no 
longer shorted out by pair production. Beyond this point, the acceleration 
grows weaker, producing fewer and lower energy Y-rays , which also cannot 
produce the particle-photon cascades necessary for the radio emission. 

The fact that there are few radio pulsars (with measured values of P) which 
have characteristic ages greater than ~ 3x10 7 years implies that the radio 
emission also turns off around this value of T and may be the result of a 
Y-ray turnoff. 

These results on pulsar Y-ray efficiency roughly agree with the estimate 
of Buccherl et al. (1978) that T]^ a t. They also have estimated Y-ray fluxes 
for 88 pulsars with measured values of P and listed those considered to be 
potentially observable. All of the pulsars in Table 1 do not appear in the 
list of Buccherl et al. because many did not have measured values of P at 
the time. They also chose a lower limit to the flux detectability of the 
COSB satellite of 0.1 times the Crab flux and therefore included some pulsars 
with lower predicted fluxes which do not appear in our list. 




V, DISCUSSION 


Thin paper has Investigated the dependence of pulsar y-ray spectra 

and luminosities on the pulsar period and magnetic, Meld strength in a 

model where the Y-rays aru produced by curvature radiation from primary 

particles. Data on the Crab and Vela pulsars were used to determine the 

primary particle energy and density t both of which nut constraints on 

theories of particle acceleration in the magnetosphere. Using these 

parameters and the calculated spectra, we have determined how Y-ray 

luminosity scales with the quantities P and B for all other pulsars. A 

formula has also been derived for the y-ray efficiency, which is found to 

1 fl 

increase with characteristic age as t . 


General y-ray luminosity formulae for pulsars have bean derived by 

other authors. Chong and Rudarmau (1980) have estimated, from considerations 

30 

of a large potential drop in the outer magnetosphere, that *'* 3x10 B^ 

P erg s at several hundred MeV. Ayasli and pgelman (1980) consider 
this outer gap model and calculate analytically the expected spectrum and 
luminosity. They conclude that L., « B and that the y-ray efficiency 

tie* " .. ia l n tin.aA'l P I r\ rtff n n i* h4 a 


scales as t' 


These papers deal with a specific model of particle 


acceleration and include the effects of pair production only on the gap 
potential or in the form of a sharp energy cutoff, We have concentrated 
here on a detailed treatment of the pair production attenuation and its 
effects on the shape of the spectrum, The results are basically independent 
of the acceleration mechanism, given that Y isv no more than weakly dependent 
on pulsar parameters. While the derived dependence of the luminosity on 
these parameters is generally consistent with the previous derivations, the 
differences reflect the effects of pair production in the magnetosphere. 


Future instruments aboard the proposed Gamma Kay Observatory (GRO) 

should be capable of detecting many wore Y-ray pulsars. The number of 

observable sources in a disk population with surface density, N Q , at a 

given luminosity, L Q , is N fl (> Sy) » i h Q N 0 ^ s y * where Sy is the y-ray 

flux above 100 MeV, Assuming a mean pulsar period of 0,5 s, and magnetic 

field strength B q « IQ 12 gauss, b Q * 4xlQ 35 a" 1, With n local ,niltu,l ‘ m,l ‘ facQ 
density of N •“ 500 kpe“ 2 (Manchester 1979), N (> S v ) ®* (lxlQ" S /S v ) photons 

ml g ® I I 

a cm ", The number of observable sources which are pulsars should be off 
order 10 with the limiting sensitivity of CGSB and 100 with GRO, The number 
which can be detected as pulsed y-ray sources , however, is somewhat less, 
as only a fraction of those will have measured radio periods and period 
derivatives. 

i 
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FIGURE CAPTIONS 

12 

Figure 1 — Calculated spectra above 100 MeV for a magnetic field of 10 

7 

gauss, Initial primary particle energy of 1 x 10 and different 
values of the pulsar period In seconds, 

12 

Figure 2 — Calculated spectra above 100 MeV for a magnetic field of 10 gauss, 
initial primary particle energy of 3 x 10^ and different values of 
the pulsar period In seconds, 

12 

Figure 3 — Calculated spectra above 100 MeV for a magnetic field of 5 x 10 
gauss, initial primary particle energy of 3 x 10' and different 
values of the pulsar period in seconds. 

Figure 4 — Fits of calculated spectra to observed spectra of the Crab and 
Vela pulsars. Data points are from Lichti et al, 1980, 

• • 

Figure 5 — Gamma ray efficiency times the quantity (P cra jj/ P) plotted 


against characteristic age. 
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